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FLIGHT  CONTROL  SYSTEM  DESIGN  FOR  A  COVIPIJTER  CONTROLLED  AIRCRAFT  WITH 


LIMITED  SENSOt^S 
Thomas  P,  Webb* 

Ab  s  tract 

A  complete  flight  control  system  for  a  small  computer  controlled 
aircraft  was  designed  using  only  yaw  rate,  heading,  lateral  load 
factor,  airspeed,  altitude,  and  rate  of  climb  feedback.  This 
multi-input  multi-output  control  problem  was  done  using  the  classical 
root  locus  technique  on  a  linearized  system  model.  The  performance  of 
the  flight  control  system  was  then  checked  using  a  12  deg r ee - o f - f r eedom 
nonlinear  simulation.  The  simulation  results  revealed  surprisingly  good 
performance,  considering  the  limitation  on  sensors. 

I .  Introduction 

The  Department  of  Electrical  Engineering  at  the  United  States  Air 
Force  Academy  is  attempting,  through  one  of  its  senior  design  courses, 
to  design,  build,  and  fly  a  computer  controlled  aircraft.  The 
Department  of  Aeronautics  was  asked  to  help  design  the  flight  control 
ivstern  to  be  implemented  by  the  on-board  digital  computer.  The  project 
involved  building  and  testing  a  wind  tunnel  model  of  the  aircraft  to 
determine  its  aerodynamic  characteristics,  performing  mass  tests  on  the 
a-  tudi  aircraft  to  determine  inertia  characteristics,  developing  a  12 
1  e g r ee - 0 f - f r eedom  nonlinear  aircraft  simulation  program,  and  designing 
t’le  actual  flight  control  system.  This  report  describes  only  the  last 
t  a  s  l<  . 

' 1 .  Ai riraft  Description 

Ttie  aircraft  acquired  by  the  Electrical  Engineering  Department  is 
an  off  tfie  shelf  hobby  radio  controlled  airplane  called  the  "Big  Stick" 
s-jiJ  by  Hobby  Shack  in  kit  form.  This  particular  aircraft  was  chosen 
f  1  r  Its  large  size  and  d  oi'  i  I  e  fi  an  d  1  i  n  g  (]  u  a  1  i  t  i  e  s  .  The  aircraft  is 

•  V,  ,1  I  ' )  r  ,  1  iS  AT  ,  Assistant  Prof.,  Dept,  n  f  Ae  r  <  i  n  a  n  t  i  c  s  ,  US  ALA  . 


III.  The  Control  Problem 


The  purpose  of  the  flight  control  system  is  to  make  the  aircraft 
fly  an  arbitrarily  specified  (and  not  necessarily  straight)  path  given 
continuous  information  on  the  current  state  of  the  aircraft  through  a 
limited  number  of  sensors.  The  design  of  the  flight  control  system  was 
formulated  as  a  multi-input  multi-output  feedback  control  problem.  The 
actual  par  ame  ters  to  be  controlled  were  specified  as  altitude  (h), 
airspeed  (V),  and  heading 
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Figure  2.  The  Control  Problem 

Referring  to  Figure  2,  the  control  problem  can  be  visualized  as 
one  of  driving  the  values  of  h,  V,  and  i  to  those  of  h^,  V^,  and 
respectively,  where  c  denotes  the  commanded  value. 

The  flight  control  syst em  makes  inputs  to  the  aircraft  by 
ad|  istiiig  the  settings  or  deflections  of  the  elevator  (6^),  throttle 
(  .),  aiie  vis  rudder  (s^).  These  control  settings  depend  on 

')  e  .'if  'se  a  s  :  I  r  erne  n  t  s  ,  wliich  contain  i  it  t  o  r  :r:d  t  i  o  n  atiouT  the  actual 
!  ‘  1 1  j  1  aircraft,  and  lire  coiri:i..inde'.l  vjiuvs  of  altitude,  airspeed, 
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Figure  3.  Control  Process 

The  coiitrol  syste.n  design  problem,  stated  simply,  was  to  find  an 
algorithm  to  convert  the  commands  and  sensor  measurements  into  the 
proper  control  settings  to  make  the  aircraft  fly  according  to  the 
commands . 

Due  to  cost  constraints,  the  sensors  aval  'able  fo''  the  project 
were  limited  to  a  yaw  rate  gyro,  a  lateral  accelerometer,  a  heading 
indicator,  an  altimeter,  and  an  airspeed  meter.  This  is  a  very  limited 
set  ot  measurements  considering  the  job  required;  tl.eretore,  it  was 
feared  that  design  of  a  satisfactory  control  system  migfit  prove 
impossible.  Tor  instance,  note  that  tliere  are  no  position  gyros.  This 
•■'■an'  that  the  con'ro!  system  is  required  to  r,;!!  the  aircraft  in  and 
;u  i  of  t.rns  w;  t  [i  no  feedback  a’h  a  t  s  oc*  v  e  r  as  to  tl  e  bank  ar.gle  nf  the 

d.rrraft.  -Mrnosr  all  'hree-axis  autopilo's  in  use  today  have  iiosition 

grvcs  for  va'v,  pi  '  r  .h  ,  and  roll  angie  me  a  s  u  r  erne  n  r  , 

-rr'  d  '  s  ;  g  n  o  :  "  jio  s  e  s  .  the  nii-  a  s  u  r  erne  n  t  s  i  v  a  i  1  ai'  1  e  r  c  '  i.  e  flight 
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1  1  1  I  r  a  !  I  Dan  t  at.  tor 


h 


altitude 


V  airspeed 

h  rate  of  climb  (to  be  derived  by  numerically 
differentiating  h) 

I  V .  Design  Pr  ocedu  re 

"Classical"  control  theory  methods  (LaDlace  tranforms  and  root 
loci)  were  used  in  the  design  of  this  control  system  as  opposed  to 
"modern"  optimal  control  techniques  (altnough  the  state-space  matrix 
r ep r e s en t a t  1  on  was  borrowed  from  modern  theory). 

Two  important  simplifications  were  made  and  carried  throughout  the 
entire  design  procedure.  The  first  wa  s  that  the  (tontrol  syst  em  wa  s 
cintinuous  ins'ead  of  discrete.  (itecail  that  tiie  control  system  is  to 
be  implemented  hv  a  digital  computer.)  This  assumption  is  not  too 
unreasonable  provided  the  cycle  time  of  the  computer  is  quite  a  bit 
fastei-  than  the  aircraft  response.  The  second  simplification  was-  that 
of  oerCn  r  sensors.  This  means  that  the  measurements  of  the  aircraft's 
state  pro'-  ided  to  the  control  system  are  true  and  uncorrupted  by  noise. 
This  .iicatijn  may  or  may  not  be  valid,  depending  on  the  quality 

t  trie  ."risors.  Xs  is  th,e  case  in  many  feedback  control  design 
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firs.gn  procedure  consisted  of  four  major  steps: 
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Figure  4.  Flight  Controi  System  Structure  (Basic) 
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R .  Simulation 
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1  .  Non  1  .  n  !  limits 

The  f  a i  1 0  w ! n  g  c  o  n  t  r  o 1  limits  were  used  in  the  simulation 
base  d  0  n  e  s  t  i  rimi  tea  a  i  r  r  a  1  t  1  ■.  m  its: 
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F  igure  9.  LaterahDsrectional  Root  Loci  (Blow-up) 
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Figure  8.  Lateral-Directional  Root  Loci 


The  simulation  results  for  the  selected  git  ins  are  presented  in  the 
next  section. 

V .  Results 

A.  Root  Loci 

In  multi-loop  feedback  systems,  such  as  the  one  being  dealt 
with  in  this  report,  the  gains  affect  the  system  poles  in  an 
interrelated  and  complex  fashion.  A  "shotgun"  (trial  and  error) 
approach  was  used  to  initially  find  a  neighborhood  of  gains  that 
appeared  to  give  reasonable  poles.  The  gains  were  then  varied  in  a 
more  systematic  fashion  to  refine  the  gain  selection.  Gains  were 
selected  on  the  basis  of  the  speed  and  stability  of  the  resulting 
poles.  This  was  done  separately  for  the  longitudinal  and 
1  a t e r a  1  - d 1 r ec t 1 ona 1  cases. 

1.  Longitudinal  case 

The  three  longitudinal  gains  selected  were: 
deg/ft/s,  ft/s/ft,  and  K^r.!  BHP/ft/s.  Figure  7  shows  a 

segment  of  the  longitudinal  root  loci  plot  that  indicates  how  the  poles 
are  affected  by  the  gains  around  the  selected  values.  Only  the  upper 
left  quadrant  of  the  complex  plane  is  shown  since  any  values  in  the 
right  half  plane  are  unstable  and  unacceptable  and  since  the  bottom 
naif  plane  iS  a  mirror  linage  of  the  top.  The  actual  values  of  the 
poles  plotted  are  contained  in  the  computer  printout  in  Appendix  E. 


The  root  loci  were  constructed  by  varying  the  gains  and  solving 
for  the  poles  after  each  change.  This  process  had  to  be  computerized. 
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Separate  programs  were  written  on  an  Apple  microcomputer  for  the 
longitudinal  and  lateral-directional  cases.  The  results  are  discussed 
in  the  next  section. 

D.  Nonlinear  Simulation 

A  12  degree  of  freedom  nonlinear  Big  Stick  simulation  program 
was  written  for  the  Burroughs  6900  computer  at  the  Air  Force  Academy  by 
Cadet  First  Class  Haniel  A.  Draeger.  A  hard  copy  is  included  in 
Appendix  D.  This  simulation  pr'svides  a  mu  c  h  more  accurate  mat  hema  t  i  c  a  I 
model  of  t  tie  airctaft  than  tfie  linearizel  equations  which  were  used  to 
d  e  ^  e  r  m i n  e  the  control  g  a  i  n  s  .  Ba  s i  •  a  1  1 y ,  i n  e  program  numerically 
integrates  Mie  nonlinear  aircraft  equations  of  mot' on  f  r  oin  Reference  1, 
modified  to  include  the  control  system,  and  plots  out  any  of  the  state 
variables  versus  time.  The  nonlinear  equations  do  not  decouple  into 
longitudinal  and  l a  t  e  r  a  1  ■  d  i  r  e  c  '  i  o  n  a  i  sets. 

The  simulation  was  run  to  see  how  the  control  gains,  selected 
under  the  linear  assumption,  would  actually  pe'^form  in  the  real, 
nonlinear  world.  The  effects  of  such  elemsMits  as  control  deflection 
limits  and  changes  in  air  density  could  be  observed.  The  simulation 
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pole  indicates  the  stability  of  its  associated  mode  (-  for  stable,  + 
for  unstable,  0  for  neutral)  and  tlie  imaginary  part  indicates  the 
oscillation  frequency  (aperiodic  if  real).  Complex  poles  occur  in 
conjugate  pairs.  The  magnitude  of  the  pole  (distance  from  origin) 
indicates  the  speed  of  the  associated  mode.  For  example,  in  Figure  6, 
poles  1,  2,  and  5  are  aperiodic.  Poles  1  and  5  are  stable  while  2  is 
unstable.  The  mode  associated  with  pole  5  will  die  out  faster  than  the 
one  associated  with  pole  1.  The  complex  conjugate  pairs  3  and  ^ 
represent  oscillatory  modes.  Mode  3  is  stable  and  4  is  unstable. 

IMAGINARY 

A 

X3A 

-X - X - 

5  1  2 

X3B 

Figure  6.  Poles  on  the  Complex  Plane 

It  can  be  shown  (Ref.  <* )  that  the  characteristic  equation  of  a 
system  described  by  Equation  3  is  det(\[I]  -  [A])  =  0  where  \  is  an 
arbitrary  scalar  number  and  [I]  is  an  identity  matrix.  The  solutions  of- 
the  equation  for  \  are  the  poles.  We  are  interested  in  the  poles  of 
the  closed-loop  system,  Equation  10.  These  can  be  determined  by 
solving  det(ni]  -  [A-FF]  )  r  0.  The  F  matrix,  of  course,  depends  on  the 


X  4A 

- >  REAL 

X  4B 


For  the  1  a t e r a  1  - d i r ec t i ona 1  model: 


[F]  = 


0 

1  .22K 


ny 


0 

.000634K 


ny 


J-.00624K 


ny 


1  -  .00624K 


ny 


1  -  .0062'»K 


K^K 

6  ^ 


ny 


[B']  . 


K^_Kr 


0 

0  ‘  0 

By  substituting  (4)  into  (3),  we  obtain 


(7) 


(8) 


X  [A]  X  +  [B]  '[^-[F]  X  f  [B'j  u^j. 


(9> 


Combining  terms  gives 

^  z  [A  -  BF]  X  +  [BB'  ]  (10) 

The  complete  linearized  system  model  (or  closed-loop  system), 
then,  consists  of  two  independent  equations  of  the  form  of  Equation  10 
--  one  for  longitudinal  motion  and  one  for  1  a t e r a  1 -d i r ec t i ona 1  motion. 
The  vectors  and  constant  matrices  have  been  defined  for  each  case. 

C.  Determination  of  Control  System  Gains 

The  core  of  the  feedback  control  problem  is  the  selection  of 
the  control  gains.  In  this  project,  that  means  finding  values  of  , 

K^,  Ky,  K  ,  ,  K^y,  and  Kyp  that  give  satisfactory  response  to 

aircraft  heading,  altitude,  and  airspeed  commands  in  the  presence  of 
disturbances  such  as  wind  gusts. 

The  response  of  any  linear  dynamic  system  is  characterized  by  the 
roots  of  Its  characteristic  equation  (also  called  system  poles  or 
eigenvalues).  Root  loci  were  used  in  this  project  to  select  ttie  gains. 
A  root  locus  IS  a  complex  plane  plot  showing  how  a  pole  varies  as  a 
gain  IS  changed.  By  Wciy  of  review,  the  sign  of  the  real  part  of  the 


algebraic  combinations  of  the  states.  Equation  1  shows  that  this  is 
true  for  the  measurement  ti.  This  is  also  true  for  V  is  closely 

approximated  by  u.  It  is  possible,  then,  in  both  the  longitudinal  and 
1  a t e r a  1  - d i r ec t i ona 1  cases  to  express  the  controls  as: 


u  =  -  [F]  X  +  [B']  u 


where  [F]  is  called  the  feedback  matrix  and  [B  ]  is  called  the  input 

j 

matrix.  The  vector  u  is  the  command  which  is  defined  as  [h  ,V  ]  for 

c  c  c 

j 

the  longitudinal  case  and  t for  the  1  a  t  e  r  a  1  -  d  i  r  ec  t  i  ona  1  case. 

From  Figure  5  and  Equation  1  and  2  the  feedback  and  input  matrices 
were  determined  to  be  as  follows; 

For  the  longitudinal  model: 


[F]  =  0  73.33K^ 


-73.33K/ 


[B'l  =  K,K|^ 


n  is  lateral  load  factor  which  is  equal  and  opposite  to 
n ong r av i t a t i ona 1  lateral  acceleration,  normalized  to  the  acceleration 
of  gravity,  g.  It  can  be  shown  from  Reference  1  and  Appendix  C  that  for 
the  linearized  approximation, 

n  -  1.222b  t  .000634p  -  .0285  r  -  .006246^,  which  is  a  linear, 
algebraic  combination  of  the  states.  (6  is  a  combination  of  the 
states  only,  since  there  are  no  commands  that  feed  into  the  rudder.) 


STATE  VARIABLE  <u,  ft/s) 

2.60  73.40  74.20  75.00  75.80 
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N  0.00 


8.00 


16.00  24.00  32.00  40.00 

TIME  (SECONDS) 
FLIGHT  HISTORY 


48. 


00  56.00  64.00 


Figure  14.  Straight  Climb  -  h 
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Figure  16.  Level  Acceleration  Simulation  -  u 


TIME  (SECONDS) 
FLIGHT  HISTORY 

Figure  17.  Multi-command  Simulation  -  'P 
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APPENDIX  D~  Simulation  Program 
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From  Reference  1,  the  linearized  1  t  e  r  a  !  -  d  i  r  er  t  i  on  a  i  equations  of 
motion  are  (using  the  definitmins  i  r  om  Table  C-2); 
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U^6  -  .  q  e  +  YgB  ^  Y  p  e  y^r  .  Y^  6^  >  Y^  6^ 


P  -  A, r  -  L.B  +  L^p  +  L^r  e  l.  6^  + 
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3  1 


SI 


r  -  B.p  =  N„B  -t-  N,  B  +  N^p  +  N  r  N,  6  N,  4 
■L  B  lo  p'^  r  6a  6r 
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^  I 

XX 


X  7 


Z1 


In  addition,  the  following  approximations  were  added: 

<3  =  P 

■  -  r 

> 

The  values  were  substituted  in  and  the  equations  were  manipulated 
algebraically  tc  obtain  first  order  matiix  form. 
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The  remaining  necessary  coefficients  and  derivatives,  listed  below, 
were  obtained  by  linearizing  the  drag  polar  around  the  steady-state 
condition  and  by  assuming  that  the  propeller  thrust  acts  through  the 
aircraft  center  of  gravity. 


Cn  =  .073 


0™  =  0 


C-j.  =  .073 


,259/rad 


C.  =  .376 
“-I 


=0 

mj 

a 


From  Reference  I,  the  linearized  longitudinal  equations  of  motion  are 
(using  the  definitions  from  Table  C-1): 

u  =  -GV+  X  u  +  Xy  u  +  X  oi  +  X,  5  +  Xy  . 

-  U^q  =  Z^u  +  Z^a  +  Z^a  +  Z^q  +  Z^^6e 


q=Mu+MyU+Ma+  My  a  +  M- a  +  M  q  +  M.  6e 
^  u  T  a  T  a 

u  a  e 


In  addition,  the  following  approximations  were  added: 


h  =  U^o  -  U^a 

The  values  were  substituted  in  and  the  equations  were  manipulated 
algebraically  to  obtain  first  order  matrix  form. 
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Estimated  Aerodynamic  Data  (Reference  1  and  5): 
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^‘wind  tunnel  data  not  used  because  model  had  no  wing  dihedral 
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fass  Data  (Reference  4): 
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APPENDIX  C-  Aircraft  Linearized  Equations  of  Motion 

The  methods  of  Reference  I  were  used  to  linearize  the  aircraft 
equations  of  motion  about  a  steady-state  condition  of  coordinated, 
straight  and  level  flight  at  7300  feet  and  73.33  ft/s  (50  mph).  The 
angle  of  attack  required  for  this  condition  was  7.2  degrees  relative  to 
body  axes.  Since  the  linearized  model  uses  stability  axes 
(longitudinal  axis  parallel  to  the  steady-state  relative  wind),  the 
inertia  terms  I  ,  I  ,  and  1  ,  of  Appendix  13  had  to  be  transformed 
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from  body  axes  to  stability  axes. 
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